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Abstract 

We use Monte Carlo simulations to study NiFe204 nanoparticles. Finite size and surface effects 
differentiate them from their bulk counterparts. A continuous version of the Wang-Landau algo- 
rithm is used to calculate the joint density of states g(M z , E) efficiently. From g(M z , E), we obtain 
the Bragg- Williams free energy of the particle, and other physical quantities. Lhe hysteresis is 
observed when the nanoparticles have both surface disorder and surface anisotropy. We found that 
the finite coercivity is the result of interplay between surface disorder and surface anisotropy. If 
the surface disorder is absent or the surface anisotropy is relatively weak, the nanoparticles often 
exhibit superparamagnetism. 



PACS numbers: 75.50.-y 75.60.-d 75.75.+a 



I. INTRODUCTION 



Many magnetic nanoparticles are found to be superparamagnets,— since they are not suf- 
ficiently big to contain more than one magnetic domain. On the other hand, a large fraction 
of the spins are located near the surface of the nanoparticle, and their local environments 
have lower symmetries than the bulk. These spins may prefer to align to directions differ- 
ent from the bulk ordering direction, thus reducing the total magnetization. This surface 
anisotropy might contribute significantly to the total magnetoanisotropy and the hysteresis 
of the nanoparticle,— e.g. ferrimagnetic NiFe2C>4 nanoparticles. The bulk NiFe2C>4 is a ferri- 
magnet with a Neel temperature Tjv = 838 K. The ball-milled nanoparticles have an average 
radius of 6.5nm^ and have been previously found to show a slightly open hysteresis loop at 
low temperatures.— £ 

In this paper, we use the heatbath algorithm and the Wang-Landau algorithm to calculate 
the finite-temperature properties of individual NiFe2C>4 nanoparticles. The heatbath algo- 
rithm 6,7 runs at a certain temperature and external magnetic field. However, the heatbath 
algorithm samples does not directly generate the hysteresis, because the canonical ensemble 
it samples is dominated by the ground state rather than a meta-stable state, which con- 
tributes to the hysteresis. We use a continuous version of the Wang-Landau algorithm^ to 
calculate the joint density of states of the nanoparticles, g(M z , E), where M and E are mag- 
netization and internal energy of the nanoparticle respectively. g(M z ,E) contains complete 
information of the nanoparticle at any temperature and in any magnetic field. 

II. MODEL OF NiFe 2 4 NANOPARTICLES AND COMPUTATIONAL METHODS 

Bulk NiFe2C>4 has an inverted spinel structure of lattice constant 8.34 A. The eight 
tetrahedral sites in one unit cell are occupied by Fe 3+ , while eight Fe 3+ and eight Ni 2+ are 
distributed on the sixteen octahedral sites. We assume the ions on the octahedral sites have 
Verwey order, i.e. the Fe 3+ and Ni 2+ occupy (100) planes alternatively. In this system, 
S = 5/2 for Fe 3+ , and S — 1 for Ni 2+ . The Hamiltonian for the magnetic structure of the 
bulk consists of Heisenberg indirect exchange: 
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where the summation is over nearest and next nearest pairs with nonzero exchange Jy. The 



exchange constants J™ are taken from Ref. 



and we use semiclassical approximation for 



the spins. In the ground state, Fe 3+ spins on tetrahedral sites are antiparallel to Fe 3+ spins 
on octahedral sites, while Ni 2+ spins align with Fe 3+ spins on octahedral sites. The cubic 
anisotropy {K\ = —8.7 x 10 4 erg/cm 3 at 77 K— ) is negligible for these nanoparticles. 

We model the nanoparticles by cutting off those spins outside a sphere of a given radius. 
The center of the sphere is randomly selected, which might not be a lattice point. Following 
the previous work 4 , we randomly remove a small fraction of surface spins to model surface 
disorder, and use uniaxial anisotropy for the surface spins: 

H A = -A; s ^(s J -n i ) 2 , (2) 

where the unit vector Sj points in the direction of spin Sj, and the unit vector fij is chosen to 
be parallel to fj) where Tj are positions of nearest neighbors of in the nanoparticle. 

As a result, the surface anisotropy term vanishes if local inversion symmetry exists for that 
spin. k s in Eq. (J2J) characterizes the strength of the anisotropy. A realistic model would have 
k s dependent on the local environment of the spin, but in our simulations, k s is a positive 
tunable constant. We use 0, 2.5 K, 300 K, and 500 K for k s , to observe the effect of different 
strength of the anisotropy. 

We use the heatbath algorithm^ and the improved Wang-Landau algorithm*^ in 
our simulations. The heatbath algorithm reaches equilibrium in presence of disorder and 
anisotropy much faster than the Metropolis algorithm.— The continuous version of the Wang- 
Landau algorithm^ calculates the joint density of states 



g{M z ,E) = J 6 



E 



M^r-^HMs), (3) 
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where N is the number of spins, S collectively denotes all the spins in the model, dfi(S) = 
Y\i S~ 2 5(Sf — S^)dS^dSfdS^ is the volume measure of the phase space, E and M z are the 
average energy per spin and average magnetization per spin. With g(M z , E) in Eq. (jSJ) at 
hand, one calculate the partition function, and the Bragg- Williams free energy defined as 
a function of temperature and magnetization: F(M Z ,T) = — ksT In J g(M z , E)e~@ EN dE. If 
the system exhibits a first order phase transition, F(M Z ,T) has double minima at constant 
T. The free energy difference between the ground state and the metastable state as well as 
the free energy barrier between them are obtained from F(M Z ,T). 



III. RESULTS AND DISCUSSION 



Using the heatbath algorithm, we calculated the specific heat and the Binder cumulant 
of magnetization for bulk samples of different sizes with periodic boundary condition. With 
finite size scaling, we obtain the Neel temperature for ferrimagnetic phase T/v = 869K, 
which is about 4% higher than the experimental value. We also compared the simulations 
for bulk material to those for the nanoparticles. Figure ^ shows the effect of the size and 
surface of the nanoparticle. 

The peak in the specific heat of the nanoparticle appears at a lower temperature, about 
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FIG. 1: Comparison between the bulk material and the nanoparticle. Left: the specific heat of 
nanoparticles has a smooth peak below the Neel temperature of the bulk material. Right: the 
reduced magnetization of the nanoparticle. The bulk sample is a cubic with 7 unit cells along each 
axis. 
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as M = N~ 



. The magnetization of the nanoparticle is substantially suppressed. 
Near the transition temperature, the magnetization of the nanoparticle gradually ramps up 
with decreasing temperature, while the bulk materials show a magnetization curve typical 
of the ferrimagnetic phase transition. 

Figure 121 shows g(M z ,E) for a typical nanoparticle with surface disorder and 387 spins. 
Figure 12(a) shows g(M z , E) for a large energy range: —1200 K< E < K, which illustrates 
the generic behavior of density of states of these nanoparticles. Here the magnetization M z 
is the average ^-component of the spins. At high energies, g(M z ,E) has a single maximum 
value at M z = for constant E. At low energies, this single maximum is replaced by a 
relatively flat part, i.e. dg(M z , E)/dM z « 0, which indicates a non-zero spontaneous mag- 
netization. The energy range of Fig. E^a) is sufficient for calculating temperature dependent 
properties for approximately T > 400 K. The magnetization as a function of temperature 
and external field is shown in Fig. E^b). A spontaneous magnetization develops at about 
T = 700 K, consistent with the reduced transition temperature shown in Fig. Since 
the model has continuous degrees of freedom, g(M z ,E) is logarithmically divergent near 
the ground state energy.^ Thus, the lowest energy that a simulation can efficiently sample 
is restricted by its resolution in energy. Fig. |2fc) and (d) show g(M z ,E) at low energies 
calculated with higher resolution in energy. Their energy ranges allow a temperature range 
180 K< T < 220 K when the g(M z ,E) is used to calculate temperature dependent quan- 
tities. In Fig. EKc), a double-valley structure appears between M = —0.4 and M = 0.4. 
This structure is responsible for a finite coercivity, because the magnetization is blocked by 
it when the magnetic field is slightly reversed. The double valley structure becomes more 
pronounced at lower temperatures. On the contrary, in Fig. Efd), the density of states of 
the same nanoparticle with a reduced surface anisotropy constant k s = 2.5 K, does not show 
any noticeable structure as in Fig. Efc). The energy range in (d) is higher than that in (c) 
because of an over all energy shift caused by the change in surface anisotropy. However the 
energy range of (d) corresponds to a temperature range of 40 K< T < 80 K. 

The lack of a finite coercivity in Fig. Efd) is due to the result of the rotational symmetry. 
Obviously the surface disorder preserves the rotational symmetry. Although there is a small 
k s , the surface anisotropy is too weak to destroy the rotational symmetry at temperatures 
much larger than k s . Consequently, the nanoparticle is well described by superparamag- 
netism. 
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When we assume a strong surface anisotropy, which is "larger" than the thermal energy 
kT and comparable to the exchange field on the surface, then we indeed observe various 
structures in g(M z , E), which are responsible for coercivity at finite temperatures. Figure El 
shows four examples of this behavior. We plot the Bragg- Williams free energy per spin 
F(M Z ,T)/N in Fig. El Clearly, all of these four examples appear to have a barrier of free 
energy between two states of spontaneous magnetization, also small oscillations on top of 
the barrier is shown in Fig. E£b) . 

To observe the effect of the intermediate surface anisotropy, we used k s = 50 K and 
performed the calculation with the same set of nanoparticles that we used for the weak 
and strong surface anisotropy. We found most of them do not appear to have a noticeable 
depression in the density of states or a peak the free energy even at low temperatures with 
a few exceptions. One of them is the nanoparticle used to do the calculation for Fig.|3fd), it 




FIG. 2: Density of states and magnetization of a typical nanoparticle containing 387 spins. The 
normalization constant in the density of states is arbitrarily selected. See the text for detailed 
explanations. 
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FIG. 3: Bragg- Williams free energy per spin, a variety of structures are observed in nanoparticles 
with strong surface anisotropy and disorder, at relatively high temperature. In each panel, four 
free energy curves are for temperatures 187 K, 216 K, 244 K, and 272 K respectively from top to 
bottom. 

was found that with k s = 50 K, its Bragg- Williams free energy for 35 K< T < 55 K appears 
to have the similar shape as shown in FigEfd). 

In summary, we have studied the magnetic properties of NiFe204 nanoparticles using 
two kinds of Monte Carlo simulations. We have found that only both surface disorder and a 
pretty strong surface anisotropy are required to develop the hysteresis. We have also noticed 
that a strong surface anisotropy is required because our model has a near-spherical shape. 
The surface anisotropics at different positions are likely to cancel each other leaving a very 
small overall anisotropy. In fact, we have observed that the total magneto-anisotropy energy 
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is only of order 0.1k s . Our approach is purely based on thermodynamics and presents no 
information on the spin dynamics of the nanoparticle. It is very likely that the surface spins 
form a disordered glass-like layer, which slows down the reversal of magnetization. 
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